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Abstract. Hereditary chip- firing models generalize the Abelian sand- 
pile model and the cluster firing model to an exponential family of games 
induced by covers of the vertex set. This generalization retains some de- 
sirable properties, e.g. stabilization is independent of firings chosen and 
each chip-firing equivalence class contains a unique recurrent configura- 
tion. In this paper we present an explicit bijection between the recurrent 
configurations of a hereditary chip-firing model on a graph and its span- 
ning trees. 



Keywords: Abelian sandpile model, chip-firing, Cori-Le Borgne algo- 
rithm, Dhar's burning algorithm, G-parking functions, recurrent config- 
urations, spanning trees. 



1. Introduction 

Chip-firing on graphs has been studied by several different communities 
over the past 25 years. In statistical physics it was introduced by Dhar 
|15j as an example of self organized criticality as proposed in the Bak-Tang- 
Wiesenfeld model [2) . In graph theory it was investigated by Tardos [M] and 
Bjorner, Lovasz and Shur [9] extending ideas introduced by Spencer [33] . 
Later Biggs [7] studied this game as related to algebraic potential theory 
on graphs. Most recently, Baker and Norine [3] have shown that using the 
language of chip-firing, one can derive a Riemann-Roch theorem for graphs 
analogous to the classical statement from algebraic geometry. 

In the Abelian sandpile model, vertices are restricted to fire individually. 
This is in contrast to the cluster firing model where vertices are allowed to 
fire simultaneously. A chip-firing model is a collection H of subsets of the 
vertex set, those subsets which are allowed to fire simultaneously if no vertex 
is sent into debt. If every vertex appears somewhere in H, and the family H 
has the hereditary property, we say that H is a hereditary chip-firing model. 
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From this perspective, the sandpile model is the coarsest hereditary chip- 
firing model, described by taking H to be the collection of all singleton sets 
from V{G) \ {vq}, and the cluster model is the finest hereditary chip-firing 
model, described by taking H to be the power set of V(G) \ {vo}. 

Some of the fundamental properties of the Abelian sandpile model and 
the cluster firing model extend to arbitrary hereditary chip-firing models: 
the stabilization of a configuration is independent of the firings chosen and 
each chip-firing equivalence class contains a unique recurrent configuration. 
It is well known that the number of chip-firing equivalence classes is the 
same as the number of spanning trees of a graph. It follows that the number 
of recurrent configurations in a hereditary chip-firing models is the same as 
the number of spanning trees. 

For the case of ASM and CFM, several bijections between recurrent con- 
figurations and spanning trees exist in the literature, e.g. [15] [2] [8] |12j . 
There is a simple relationship between the recurrent configurations in ASM 
and CFM which allows a bijection is one model to be "dualized" to produce 
a bijection in the other model. The recurrent configurations in CFM go by 
several names: G-parking functions, uo-reduced divisors, superstable config- 
urations. It is the aim of this paper to present an explicit bijection between 
the recurrent configurations in an arbitrary hereditary chip-firing model and 
the spanning trees of a graph. Our bijection is a modification of the Cori-Le 
Borgne algorithm |14j . 

If we order the elements of H by inclusion, we have a set of maximal 
elements Ai, . . . Af. which in turn, by the hereditary property, determine H. 
We note that these maximal elements of H need not be disjoint, i.e. hered- 
itary chip- firing models are not determined by partitions of the vertex set, 
instead we should think of them as being described by covers of V(G) \ {vq}. 
Moreover, these covers need not be irredundant, instead we ask that the 
elements of the cover be incomparable. This allows us to naturally identify 
hereditary chip-firing models with maximal antichains in the Boolean lattice. 
Calculating the number of such maximal antichains is a challenging problem 
[17] . but this quantity is easily seen to be at least exponential in n. 

2. Notation and Terminology 

We take G to be a connected undirected loop less multigraph with vertices 
labeled v , Vl ,...v n . Given X,Y C V(G), we let (X, Y) = {e G E{G) : e = 
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(x,y),x £l,f/£ Y}, and let X c denote V(G) \ X. To describe chip-firing, 
we begin with a graph G and a configuration D of chips on G. Formally, a 
configuration of chips is a function D : V[G) — > Z n+1 . For the purposes of 
this paper we will usually restrict our attention to D such that D{vi) > 
for all i ^ and D(vo) = — Ya=i D{vi) so that the sum of the values of 
D, called the the degree of D, is 0. If a vertex v in a configuration of D 
is seen to have D(v) < 0, we say that this vertex is in debt. The basic 
operation is firing whereby a vertex v sends a chip along each of its edges to 
its neighbors and loses deg(v) chips in the process so that the total number 
of chips is conserved. We designate vq to be the sink vertex and say that it 
cannot fire. This ensures that we cannot continue firing vertices indefinitely. 
The adjacency matrix A of a graph isan+lxn+1 matrix with entries 
■A-ij = # of edges between Vi and Vj. Taking D to be the diagonal matrix 
with Di t i = degree of the Laplacian of a graph is defined as the difference 
D — A. 

For S C V(G), we take xs to be the characteristic vector of S. As an 
abuse of notation we denote X{vi} by Xi- From a linear algebraic perspective, 
viewing a configuration D as a vector, if a vertex Vi fires then D is replaced 
by D — Qxi, and more generally if a set S fires we obtain D — Qxs- We say 
that two configurations D and D' are equivalent if there exists some sequence 
of firings which brings D to D' (possibly including firings by vq and passing 
through intermediate configurations which are negative at vertices other than 
vq). Two configurations are seen to be equivalent if their difference is in 
the integral span of the columns of the Laplacian. We call a collection of 
configurations which are equivalent, a chip-firing equivalence class. 

The ASM (Abelian sandpile model) is defined by placing the additional 
restriction that vertices may only fire one at a time, whereas in the CFM 
(cluster firing model), vertices are allowed to fire simultaneously. We fix a 
collection H of subsets of V(G) \ {t>o}, those sets which are allowed to fire 
simultaneously if no vertex is sent into debt, and call this collection a chip- 
firing model. If each vertex vi with i / appears somewhere in H, we say 
that H covers G. If H covers G and H is hereditary, i.e. for every A G H 
and B C A, we have that B G H, we say that H is a hereditary chip-firing 
model. 

Let H be a hereditary chip-firing model on a graph G. If a configuration 
of chips D has no set of vertices M £ H which can fire without some v G M 
being sent into debt, we say that D is stable. The process of firing sets from 
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H until a configuration becomes stable is called stabilization. We say that 
a set M 6 -ff, is ready in D if this set can fire without sending any vertex 
into debt, and call a vertex v active in a configuration Z) if there exists some 
M C V(G) \ {vo} with v € M which is ready. Suppose t> € V(G) is active in 
a configuration D. There may very well be several different maximal ready 
sets which contain v, and these different maximal ready sets might cause 
v to lose different numbers of chips if they were to fire. Therefore, we let 
m(v,D) denote the minimum amount that an active vertex v can lose by 
firing a maximal ready set in D which contains v. 

Lemma 1 states that the stabilization of a configuration in a hereditary 
chip-firing model is well defined, so we denote the stable configuration ob- 
tained from D by stabilization as D°. A configuration D is said to be reach- 
able from another configuration D' if there exists a way of adding chips to D' 
and then firing ready sets to reach D. Because of our convention that the de- 
gree of D be zero, we are actually adding configurations of the form Xi — Xo-, 
i.e. subtracting from vq exactly as many chips as we add to other vertices. 
This operation is referred to in [3] as the Abel-Jacobi map. A configuration D 
is globally reachable if it is reachable from every other configuration. Finally, 
we call D recurrent, if it is both stable and globally reachable. The original 
motivation for this terminology comes from the observation that if we con- 
tinue adding chips and stabilizing, the configurations we will see infinitely 
many times are the recurrent ones. The recurrent configurations in CFM 
(G-parking functions) are precisely the stable configurations, so there is no 
need for a discussion of global reachability. We say that a configuration u is 
critical if it is stable and (D — Qxo)° = D. As with the ASM, a configuration 
is recurrent if and only if it is critical. This statement is trivially true for 
the CFM. 



3. Preliminary Results 

In this section we present the basic results of hereditary chip-firing mod- 
els. Hereditary chip-firing models as well as the results of this section were 
discovered independently of the author by Paoletti [27] [2H], and Caracciolo, 
Paoletti and Sportiello|llj. They observe that stabilization in a chip-firing 
model H is independent of firings if and only if H is closed under subtrac- 
tion, i.e. for all A, B € H , we have A \ B € H . They then restrict to the 
case where for each v £ V(G), {v} € H. It is easy to see that a family of 
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subsets of [n] is closed under subtraction and contains all singletons if and 
only if it is hereditary and covers [n] . 

Lemma 1. Given a fixed hereditary chip-firing model H on a graph G, and 
a chip-firing configuration D on G, the stabilization of D is independent of 
the firings chosen. 

Proof. First, we observe that if M,N C V(G) \ {vq}, M is ready and N 
fires first, then M \ N is ready. This is because if we fire N and then fire 
M \ N, a vertex v 6 M loses at most as many chips as if M had fired 
alone. More generally, if M is ready and a multi set N fires, i.e. we fire 
vertices in N a number of times corresponding their multiplicity in N, then 
M \N is ready. Let M 1 ,...,M S € H and iV"i, . . . , N t G H correspond to 
sequences of sets which are fired in two different stabilizations of D. Let 
x M q = J2i=iXM t and X Br = Y7i=iXB r - Suppose that D - QX Ms and 
D — QX^ t are not equal, i.e. the two stabilizations of D are different. We 
note that this can occur if and only if Xm s ^ X^ t , as v$ does not fire and 
the kernel of the Laplacian is generated by the all one's vector. It follows 
that, without loss of generality, there exists some I maximum such that 
Xmi < Xpf t and Xm 1+1 ^ -XjVf By construction Mz +1 is ready for D — QXm v 
Now let xp = X]y t — Xmi be the characteristic vector corresponding to the 
multi set P . By the first observation, M/ + i \ P is nonempty and ready for 
D — QXu x — Qxp = D — QX]y t , but this contradicts the fact that D — QX^ t 
is stable. 

□ 

Theorem 1. Given a fixed hereditary chip- firing model H on a graph G, 
there exists a unique recurrent configuration v in each chip- firing equivalence 
class. 

Proof. We begin by observing that every chip-firing equivalence class con- 
tains at least one recurrent configuration. In a stable configuration, each 
vertex v has at most deg(u) — 1 chips. Therefore, if we can show that each 
equivalence class contains a configuration with more than deg(u) chips at 
each vertex v, it would follow that this configuration is globally reachable 
and hence its stabilization is recurrent. The technique which we now apply 
also appears in [3]. Partition the vertices according to their distance from 
vq. Let d be the maximum distance of a vertex from vq. Begin by firing 
all of the vertices of distance at most d — 1 from vq . This has the effect of 
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sending money to the vertices of distance d. Repeat until each such vertex 
v has at least deg(t>) chips. Now fire all of the vertices of distance at most 
d — 2 from vq until the vertices of distance d — 1 have at least their degree 
number of chips. Working backwards in this way towards uq, we obtain the 
desired configuration. 

We now show that there is at most one recurrent configuration in each 
equivalence class. This proof is identical to the argument presented in [20J - 
it is a little surprising that their argument requires no modification. First, 
we would like to show that there exists a configuration e with e(vi) > 
for all i ^ 0, such that when we add e to a recurrent configuration v and 
stabilize, we obtain v. Let D be a configuration such that D(vi) > deg(wj) 
for all i 7^ 0. We will take e = D — D°. Because v is recurrent, it is globally 
reachable, hence there exists some configuration £ such that [D + £)° = v. 
We are interested in computing 7 = (D + £ + e)°. Because stabilization is 
independent of firings chosen, we can stabilize 7 by first stabilizing D + £, 
i.e. 7 = {{D + C)° + e)° = {y + e)°. On the other hand, this is also equal to 
(D° + C + e)° = (D° + C + D - D°)° = (C + D)° = v. 

Assume that there are two different equivalent recurrent configurations v 
and v' such that v ~ u'. By definition, there exists some / E Z n+1 such that 
v — v' = Qf, moreover we can take / to be such that f(vo) = because 
the all ones vector is in the kernel of Q. Let / + ,/ _ € Z n+1 be such that 
/ + > 0, f~ < 0, and / + + /"" = /• Therefore, there is some configuration 
D such that D = v — Qf + = v' — Q(—f~). Note that because v and v 1 
are stable, it follows that D may have vertices which are in debt. For any 
k € N, v + ke and v' + ke will stabilize to v and v' respectively, as was 
shown above. On the other hand, if we take k to be sufficiently large, we can 
perform firings defined by / + and — /~ (by individual vertices for example) 
to v + ke and v' + ke respectively to obtain the configuration D + ke. But 
now we arrive at the contradiction that D + ke should stabilize to both v 
and v 1 . 

□ 

Lemma 2. Given a fixed hereditary chip-firing model on a graph G, a chip- 
firing configuration v on G is recurrent if and only if it is critical. 

Proof. Suppose first that v is recurrent, but not critical, that is {y — Qxo)° = 
D^u. Then (u + ke - Q X o)° = (0 + ke)° - Q Xo )° = (y - Q X o)° = D. 
Because e(vi) > for all i / 0, we can take k sufficiently large so that 
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(y + ke — Qxo){vi) > deg(fj) for all z 7^ and it follows that D is recurrent, 
a contradiction. Conversely, suppose that D is not recurrent, but that D is 
critical, then (D — kQxo) = D for all k € N. If we take k to be sufficiently 
large, then we can perform firings as in the beginning of Theorem 1 to spread 
the chips around in the graph and reach a configuration which has at least 
degree number of chips at each vertex. It follows that D is globally reachable, 
hence recurrent, a contradiction. 

□ 

Lemma 3. The number of chip-firing equivalence classes on a graph G is 
the same as the number of spanning trees of G. 

Proof. Let Q denote the matrix obtained from Q by deleting the first row 
and column, i.e. the row and column corresponding to Vq. This matrix, 
called the reduced Laplacian of a graph is known to have full rank as G is 
connected, and by the matrix tree theorem det(Q) is equal to the number of 
spanning trees of G [21j . By ignoring the values of vq in our configuration, we 
see that the number of different chip-firing equivalence classes is the number 
of cosets for the image of Q and this index is given by det(Q). □ 



4. BlJECTION 

This algorithm is a modification of the Cori-Le Borgne algorithm |14| 
as presented in [2]. Their algorithm can be viewed as a variant of Dhar's 
burning algorithm [16]. We will need to call Dhar's burning algorithm as 
a subroutine, so we first begin by describing this method, and do so in the 
context of the cluster firing model where the author believes it is more nat- 
urally understood. One might argue that the brilliance of Dhar's algorithm 
is that its discovery occurred in the context of the Abelian sandpile model, 
where its application is less obvious. 

Given a recurrent configuration v for the sandpile model K + — v = v is a 
recurrent configuration in the cluster firing model, where K + (v) = deg{v) — 1 
for all V{G) \ {vq}. The interested reader can prove this fact for themselves 
using Lemma 3 or look to [3| for an alternate proof. This allows a bijection 
for one model to be "dualized" to produce a bijection for the other model. 
The bijection presented here is the first bijection which the author is aware 
of that applies directly to both models without exploiting this duality. 
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As was mentioned in the introduction, the recurrent configurations in the 
cluster firing model are precisely the stable configurations, therefore, to check 
that a configuration v is recurrent, we need only check that there exists no 
set A C V(G) \ vq which can fire without sending a vertex into debt. A 
priori we would need to check an exponential number of sets to be sure that 
v was reduced, but Dhar's observation is that it's sufficient to check only 
n such sets. Assume that v is stable and begin by firing A\ = V(G) \ vq. 
By assumption, there exists at least one vertex v which is sent into debt. 
Remove v from A\ and continue firing sets in v and removing vertices sent 
into debt until reaching the empty set. 

Here is why this works: suppose that B £ V(G) \ vq is ready in v, but 
that we have a collection A±, . . . , A n of sets which were obtained from a run 
of Dhar's algorithm. There exists i maximum such that B C A{. It follows 
that Ai^i = Ai \ v, with v G B, where v was sent into debt by Ai, but if 
we fire Ai\B, v may only gain chips, and v is supposedly able to fire in B 
without being sent into debt. Firing Ai\B and then B is the same as firing 
Ai, contradicting the fact that v was sent into debt by Ai. 

Dhar's burning algorithm earns its name from the following alternate de- 
scription: Place D(v) firefighters at each vertex, and start a fire at vq. The 
fire spreads through the graph along the edges, but is prevented from pass- 
ing through vertices by the firefighters located there. When the number of 
edges burned incident to a vertex is greater than the number of firefight- 
ers present, the firefighters are overpowered and the fire burns through the 
vertex. A configuration is stable in the cluster firing model if and only if 
the fire burns through the entire graph. Dhar noticed that by burning in a 
systematic way, this algorithm can be turned into a bijection between the 
recurrent configurations and the spanning trees. 

In the Cori-Le Borgne algorithm, the edges are burned in a different or- 
der to produce an "activity preserving" bijection. To describe the Cori- 
Le Borgne algorithm, we begin with an arbitrary ordering of the edges 
ei,e2, • • • ,e m <G E(G). The setup is the same as with Dhar, except that 
we burn one edge at a time, always taking the edge with the smallest la- 
bel connecting the burnt vertices to the non burnt vertices. When an edge 
picked/burned causes the firefighters at a vertex to be overpowered and the 
vertex to be burnt, we mark this edge. It is clear that if the fire burns 
through the graph, these marked edges form a spanning tree . Converesly, if 
we start with a tree and begin burning the edges of our graph one at a time, 
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the edges of the tree tell us when we should burn a vertex, hence how many 
firefighters (chips) a vertex should have. This shows that the algorithm pro- 
duces a bijection between the recurrent configurations and spanning trees. 

We now give a third characterization of recurrent configurations. This 
definition is the the one which will be used in our bijection. 

Lemma 4. A configuration v is critical if and only if any maximal sequence 
of firings by active vertices brings v — Qxo back to v. 

Proof. Here, we are allowing active vertices to fire even though this may 
cause them to go into debt. If a configuration v is critical, it is clear that 
we can continue firing active vertices in the ready sets and eventually return 
to v. Conversely, suppose that there exists some firing of individual active 
vertices which brings v — Qxo back to v, but that v is not critical. If this is 
the case, there must be some vertex v G V(G) \ {vo} which was never fired 
in the stabilization of v — Qxo- We might take v to be the first such vertex, 
but observe that this situation may only occur if a vertex of the same type 
has already been fired causing v to become active, a contradiction. □ 

We let AT denote the set of recurrent configurations from a hereditary 
chip-firing model on a graph G, and T denote the set of spanning trees of 
G. Recall m(v, D) is minimum amount that v can lose by firing a ready 
set in D which contains v. We now explain our bijection between recurrent 
configurations in a fixed hereditary chip-firing model if on a graph G and 
the spanning trees of G. Here is the algorithm a for taking a recurrent 
configuration D and producing a spanning tree T = a{D): 

• Begin with X = v , R = 0, T = 0. 

• Take e\ G (X, X c ) with i minimum such that e« ^ R. 
Let ej = (u, v) with u G X and v G X c . 

• If D(v) < m(v,D - Qxx) - \{e e R : e = (w,v)} U {ej}|, set 
R := RU {a}. 

• If D(y) = m(v,D — Qxx) — |{e G R : e = (w, v)} U {e,}|, we set 
T:=TU {e;}, R := 0, and X := X U {v}. 
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We now describe our algorithm 7 for taking a tree T and producing a 
recurrent configuration, j(T). This process has two parts. First we use T 
to construct a total order on the vertices: 

7 :T->M 
Part f: 

• Begin with Xq = {vq} = {wq}, and R = 0. 

• Take ei G (X, X c ) with e% ^ R and z minimum. 

• If a £T, set i?:=i?U{ei} 

• If ei G T with e$ = («, u) and v G A" c , set v = Wi, X{ = X^\ U {wi}, 
and i? := 0. 

This process terminates when we have observed every edge of T, at which 
point we have acquired a total order w < w 1 < • • • < w n on the vertices. 

The idea is to now run the algorithm "backwards" using the total order to 
reconstruct D: 

7 : r -> M 
Part2: 

• Begin with Yq = {wq, . . . , w n -{\ and R = 0. 

• Assume that at the ith iteration of step 1, the values of D(wj) for 
j > n — i + 1 have already been calculated. 

• Step 1: Fire Y{ = {wq, . . . , w n -i-i}, i.e. compute D — Qxy i - 

• Find the maximal ready sets in V(G) \ Yi for D — QxYi by apply- 
ing Dhar's burning algorithm to Aj n (V(G) \ Yi) for each maximal 
element Aj G H. 

• Compute m(w n ^i, D - QxyJ- 

• Step 2: Pick e& with minimum so that e k ^ R. Let = (it, v) G 

• If ejfc ^ T, set R := RU {e^}, and go to step 2. 

• If efc G T, we know that e& = (u,w n -i). 

Set D(w n -i) = m(w n -i, D - Qxy,) ~ \{e e R : e = (v,w n -i)}U{e k }\, 
R := 0, and Yi := l^+i, then go to step 1. 
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Theorem 2. The operations, a and 7 are inverse to each other and induce 
a bijection between the recurrent configurations of a hereditary chip-firing 
model H on a graph G and the spanning trees of G. 

Proof. First we claim that 7 o a is the identity map on the recurrent config- 
urations. Let D be recurrent and o~(D) = T a spanning tree. Observe that 
the total order produced on the vertices of G during the run of 7 on T is 
the same as the order in which the vertices are processed during a run on 
D. Given this total order on the vertices, the algorithm 7 is designed so as 
to produce the configuration D such that a{D) = T. It follows that a is 
injective, and by Lemma 3, a is an injective map between two sets with the 
same cardinality. It follows that a is a bijection with explicit inverse 7. □ 
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